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Abstract 



> 

The formulation of an exact functional renormalization group equation for Quantum Ein- 
. stein Gravity necessitates that the underlying effective average action depends on two 
metrics, a dynamical metric giving the vacuum expectation value of the quantum field, 



and a background metric supplying the coarse graining scale. The central requirement 
of "background independence" is met by leaving the background metric completely arbi- 
trary. This bimetric structure entails that the effective average action may contain three 
classes of interactions: those built from the dynamical metric only, terms which are purely 
background, and those involving a mixture of both metrics. This work initiates the first 
study of the full-fledged gravitational RG flow, which explicitly accounts for this bimet- 
ric structure, by considering an ansatz for the effective average action which includes all 
three classes of interactions. It is shown that the non-trivial gravitational RG fixed point 
central to the Asymptotic Safety program persists upon disentangling the dynamical and 
background terms. Moreover, upon including the mixed terms, a second non-trivial fixed 
point emerges, which may control the theory's IR behavior. 



1 Introduction 



Background independence constitutes one of the central guiding principles in the quest 
for a viable quantum theory of gravity. This requirement is central in loop quantum 
gravity [1, 2, 3] and also implemented in lattice approaches towards quantum gravity [4]- 
[7] . Loosely speaking, it implies that the spacetime structure realized in Nature should not 
be part of the theory's definition, but rather emerge from a dynamical principle. This strict 
background invariance, referring to no background structure whatsoever, is, however, very 
hard to implement. In particular, without an ab initio metric the notions of causality and 
equal time commutation relations are not defined, so that the usual quantization procedures 
underlying ordinary quantum field theories cannot be applied straightforwardly. 

A milder, but nevertheless equally admissible road towards a viable quantum grav- 
ity theory is the requirement of "background covariance". This allows to introduce a 
background metric as an auxiliary tool, as long as none of the theory's basic rules and as- 
sumptions, calculational methods, and predictions, depend on this special metric. In other 
words, all metrics of physical relevance are obtained from the dynamics of the theory. This 
is the viewpoint adopted in many continuum field theory approaches to quantum gravity, 
in particular by the functional renormalization group approach initiated in [8]. 

The latter is based on a functional renormalization group equation (FRGE) which 
encodes a kind of (continuous) Wilsonian RG flow on the space of diffeomorphism invariant 
action functionals. These actions naturally depend on the expectation value g^ u = (7 M „) of 
the quantum metric 7^. In addition, the coarse graining operation requires a background 
structure, which can be used to define volumes over which the quantum fluctuations are 
averaged. This structure is conveniently provided by the background field method [9] which 
also ensures the background covariance of the approach. Here the quantum metric is split 
according to 

lnu = 9^u + (1.1) 

where is a fixed, but unspecified, background metric and h^ u are the quantum fluc- 
tuations around this background which are not necessarily small. This allows the formal 
construction of the gauge-fixed (Euclidean) gravitational path integral 

J VhVC^VC, ex P {-S[g + h] - S s % g] - S* host [h, C, C; g] - A k S[h, C, C;g}}. (1.2) 

Here S[g + h] is a generic action, which depends on 7^ only, while the background gauge 
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fixing S gi [h;g] and ghost contribution S ghost [h,C,C;g] contain g^ u and in such a way 
that they do not combine into a full 7^. They have an "extra i^-dependence" and are 
not invariant under split-symmetry 5h^ v = e^, Sg^ v = —e^ u which is respected by the 
combination (1.1). 

The key ingredient in the construction of the FRGE is the coarse graining term A k S[h, C,C;g]. 
It is quadratic in the fluctuation field h, f d d Xyfg h fJiU 7Z k lupcr (— D 2 )h f „ , plus a similar term 
for the ghosts. The kernel IZf^^ip 2 ) provides a ^-dependent mass term which separates the 
fluctuations into high momentum modes p 2 ^> k 2 and low momentum modes p 2 <C k 2 with 
respect to the scale set by the covariant Laplacian of the background metric. The profile of 
lZ l k lI/pa (p 2 ) ensures that the high momentum modes are integrated out unsuppressed while 
the contribution of the low momentum modes to the path integral is suppressed by the 
/c-dependent mass term. Varying k then naturally realizes Wilson's idea of coarse graining 
by integrating out the quantum fluctuations shell by shell. 

Taking the formal ^-derivative, eq. (1.2) provides the starting point for the construction 
of the functional renormalization group equation for the effective average action T k [10, 
11]. (See [12] for reviews.) For gravity this flow equation takes the form [8] 



Here t = \og(k/ko), STr is a functional supertrace which includes a minus sign for the 

ghosts £ = (C),£ = (C), TZk is the matrix valued (in field space) IR cutoff introduced 

(2) 

above, and T k is the second variation of Tk with respect to the fluctuation fields. Notably, 



The explicit dependence on the two metrics is essential for being able to write down the 

(2) 

exact flow equation (1.3), as the Hessian T k is the variation of T k with respect to the 
fluctuation fields at fixed g^ u . In this sense, Tk and its flow is of an intrinsically bimetric 
nature. In particular, the construction of T k involves the terms A k S and S gi + S ghost where 



terms therefore provide a source for the extra background field dependence of IV To stress 
this point, it may be illustrative to write 




(1.3) 



T k [h, £, £; g] depends on two metrics, the background metric g^ v and the expectation value 




the (^-dependence does not combine with into the full averaged metric g^. These 



F k [g,g,£„£\ = F k {h = g - g,£,£;g], 



(1.5) 
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where Tk now depends on two full fledged metrics, g^ and g^ v . 

One way to extract physics information from the FRGE is by applying perturbation 
theory [13, 14, 15]. The main virtue of the flow equation is, however, that its use is 
not limited to perturbation theory and can also be employed to obtain non-perturbative 
information. Here the most common approximation scheme consists of truncating the 
space of functionals Tk to a finite-dimensional subspace and projecting the flow equation 
onto this subspace. Studying the gravitational RG flow, within these truncations, the most 
exciting result obtained to date is a substantial body of evidence [8], [13]- [50] in support of 
Weinberg's asymptotic safety scenario for gravity [16, 17, 18]. All truncations of the FRGE 
have displayed a non-Gaussian fixed point (NGFP) of the gravitational RG flow and there 
is also mounting evidence [13, 14, 32, 33] that its number of relevant couplings is actually 
finite. This fixed point may thus provide a fundamental and predictive UV completion 
of gravity within Wilson's generalized framework of renormalization. (See [51]- [54] for 
cosmological applications of this framework.) 

While already impressive, a serious caveat in this body of evidence is that all compu- 
tations carried out to date are essentially "single-metric" and do not properly reflect the 
bimetric nature of the flow equation. Typically the ansatz made for Tk falls into the class 



where T k [g] are interaction monomials built from the expectation value metric only. The 
split-symmetry violating interactions built from both h, g are encoded in T k [h; g] which, by 
construction, vanishes for h — 0. Finally, S gi and S ghost are taken as the classical gauge- 
fixing and ghost terms. The single-metric computations then proceed by taking the second 
variation of T k with respect to the fluctuation fields and setting g = g afterwards. This 
suffices to extract the running of the coupling constants contained in r\[g]. 

The potentially problematic feature of these computations is that the /3-functions en- 
coding the running of the coupling constants multiplying interactions built from the "gen- 
uine" metric g^ v are tainted by contributions originating from pure background terms. A 
single-metric truncation does not distinguish between the running of, say, 



r k [h, £, 6 g] = f k [g] + r k [h; g\ + s gf + 



(1.6) 




and 




(1.7) 
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It determines the running of a linear combination of the couplings G%, and their back- 
ground analogs G®, A® only. 

The bimetric nature of the gravitational average action has been appreciated only very 
recently, by carrying out preliminary studies in conformally reduced gravity [43] , and study- 
ing the bimetric terms induced by quantum effects in the matter sector [44] . Currently there 
are no results on full-fledged gravity available. There are, however, several good reasons 
why disentangling between the g^ v and g^ v contributions is of central importance. Firstly, 
the single- metric computations do not account for the bimetric nature of the gauge-, ghost- 
, and cutoff-terms which inject an extra ^-dependence into the path integral. They are 
sources of split-symmetry breaking action monomials, which will inevitably be "switched 
on" along the RG flow, leading to new interactions which are either constructed from the 
background metric only or a mixture of background and expectation value metric. The 
preliminary results obtained in [43] and [44] suggest that disentangling these interactions 
may lead to a significant alteration of the results obtained in the single-metric case. In 
particular, separating the g and g pieces in Y k may destroy the NGFP underlying Wein- 
berg's Asymptotic Safety idea. Secondly, identifying g = g does not probe the direction 
of the IR cutoff AkS in theory space which may give rise to a important contribution to 
the RG flow in the UV. Thirdly, employing the background field method, the counterterms 
found in perturbation theory are constructed solely from the background fields [55], so that 
isolating their effect requires careful distinction between the g^ v and g^ v field monomials. 

Based on this motivation, our work initiates the first study of the full-fledged gravita- 
tional RG flow in a fully bimetric setting. Concretely, we study the RG flow in the bimetric 
Einstein-Hilbert truncation which distinguishes the four action monomials in eq. (1.7), i.e. 
\J~9i V9, \[9 Ri \/9 R-i respectively with R = R{g). As our central result, we show that 
the non-Gaussian fixed point known from the single-metric computations also appears in 
the bimetric case. Subsequently, we supplement the bimetric Einstein-Hilbert truncation 
by a prototypical one-parameter family of mixed action monomials built from both the 
"genuine" and background metric. Including the extra interactions induces a split of the 
known NGFP in a UV and an IR fixed point, which can be connected by a complete RG 
trajectory. 

The rest of the paper is organized as follows. In Section 2 we describe the details of 
the setup and state our main new result: the /3-functions of the double-Einstein-Hilbert 
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truncation in four dimensions. The properties of these /3-functions are analyzed in Section 
3 and we discuss our findings in Section 4. A brief summary of the heat-kernel techniques 
employed in the paper and the rather lengthy /3-functions for the bimetric Einstein-Hilbert 
truncation valid for any spacetime dimension d are relegated to the Appendices A and B, 
respectively. 

2 ^-functions of the double-Einstein-Hilbert trunca- 
tion 

In this section we derive the /3-functions of the double-Einstein-Hilbert truncation. Besides 
the Einstein-Hilbert action constructed from g^ v , known from previous single-metric trun- 
cations, the corresponding truncation ansatz also encompasses a Einstein-Hilbert action 
constructed from the background metric g^ v and a simple class of interaction monomials 
including both g^ u and g^ u . This prototypical setup accounts for the bimetric nature of 
the FRGE (1.3), for the first time disentangling the quantum gravity effects in and g^ u 
in a full gravity computation. 1 

2.1 The truncation ansatz 

Our ansatz for the double-Einstein-Hilbert truncation takes the form 

F k [g, C, C, g] = r- ct [g, g] + S* [h- g] + S^ host [g, C, C, g] , (2.1) 

where T™ et [g,g] is the metric part of the effective action (built from both g and g) which 
we supplement by the classical gauge-fixing and ghost action S gf and S shost , respectively. 
Explicitly, we consider the following one-parameter class of gravitational actions 

Here the unbared (bared) quantities are constructed from the expectation value metric g a p 
(background metric g a fi)- Furthermore, Gt and A& denote the Newton's constants and 

1 For a related analysis in the framework of conformally reduced and matter induced gravity, see ref. 
[43] and [44], respectively. 



J d d x^g [R - 2A*] 
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cosmological constants, with the superscript A and B indicating that the corresponding 
interaction term is constructed from g a p and g a p respectively. The form of the bimetric term 
appearing in the last line is motivated by the structure of the flow equation encountered in 
[44], where it is precisely the ratio y/g/y/g that naturally appears on its right-hand- side. 
In the following, we will consider integer exponents n > 2 only. The terms n — 0, n — 1 
give rise to the monomials multiplying A^ or A®, respectively, and are already included in 
the Einstein-Hilbert actions constructed from the "genuine" and background metric. 
In the sequel, we will work with the geometric gauge-fixing, setting 

Sgi = hJ ^^F^K , F, = D»h, u - \D^h , (2.3) 

and subsequently taking the Landau gauge limit a — > 0. As it was shown in [13, 32], this 
gauge choice is perfectly adapted to the transverse-traceless decomposition utilized in Sec- 
tion 2.3 below, where it leads to significant simplifications. The ghost action exponentiating 
the resulting Faddeev-Popov determinant takes the form 

^ghost = J d d xy/^C ll M^ u C u , (2.4) 

with 

M , u = gWg<*D x { 9pu D a + g av D p ) - I g^r X D x g au D p . (2.5) 
This completes the specification of our truncation ansatz. 

2.2 The conformal projection technique 

Our next task is to project the gravitational RG flow onto the subspace spanned by (2.1), 
so that we can compute the /3-functions for the /c-dependent couplings contained in T™ et . 
Obviously, this cannot be achieved by evaluating the flow equation setting g a p = g a p, which 
underlies the single-metric computations. Instead, we will resort to the conformal projec- 
tion technique introduced in [44] which identifies g^ u and g^ u up to a constant conformal 
factor: 

g pu = {1 + e) u g^ , v = ^ , f = 1 + e . (2.6) 
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Substituting this identification into (2.2) and performing a double-expansion in e and the 
background curvature R yields 



-pmct 
1 k 



[9,9] 



gF 



+ 



(d-2) 



Afe 



(1-n) 



gF 



+ 



gF 



(l-n)(2-rarf) M fc 
2 G*j 



(2.7) 



4r + — 

G? Gf 3 



+ GF e 



+ 



where the dots indicate higher powers in the e, R expansion. Plugging this expansion into 
(1.3), the left-hand-side of the equation indicates that the running of the coupling constants 
contained in the ansatz (2.2) is captured by the coefficients 



(2- 



RP . , e 1 , e 2 , 

of this double expansion. Thus, by extracting the corresponding contributions from the 
right-hand-side of the flow equation, we are able to disentangle the running of G A , A A and 
Gf , together with M k . 

At this stage, we feel obliged to add the following word of caution. While the conformal 
projection technique employed here is capable of distinguishing between the running cou- 
pling constants associated with monomials built from g^ u , the background metric g^ u or a 
mixture of the two, it has only limited power for resolving different tensorial structures. 
As an illustrative example, we consider the following three mixed terms in d = 4: 



/ 



d x (y/gy/g) 



d 4 x^g(g, v g^), 



d 4 xVd(r u g^). 



(2.9) 



Under the conformal identification (2.6) all three invariants are projected onto the same 
structure, £ 2 J d^Xy/g, and are thus indistinguishable. 2 Resolving this ambiguity will re- 
quire a much more sophisticated computational technique, like the /i Q( g-expansion advo- 
cated in [44]. Owed to the increased technical complexity of the /^-expansion, however, 
we will refrain from resolving this ambiguity and resort to the conformal projection scheme 
in the sequel. In any case, we expect that the latter is sufficiently elaborate to give some 
first insights into the properties of the gravitational RG flow taking the bimetric nature of 
the flow equation into account. 



2 The situation is completely analogous to the projection of the three i? 2 -couplings on a spherical 
background [13, 22, 23, 24, 33], which also determines the /3-functions for one particular linear combination 
of the three couplings only. 
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2.3 Hessian cutoff implementation, and the flow equation 

In the next step, it is convenient to first compute the quadratic forms arising at second 
order in the /i^-expansion of T™ et . These forms will considerably simplify the computation 
of the Hessian later on. We start by constructing the Taylor series of r™ et around the 
background (2.6): 

rn^] = rpfag] + ^T'bM^ + \5 2 vT[g,g]\ g =^- g + .... (2.10) 

To simplify the notation it is useful to abbreviate the interaction monomials in T™ et by 

h = Jd d x^R, l2 = j^ x ^{^jj) ■ (2 ' n) 

The interaction term multiplying A^ and are special cases of I 2 , corresponding to 
n = 0,1. 

Expanding the curvature invariants up to second order in h^ u and setting g^ u = l 2u g^ v 
afterwards, the second variations become 

5% = J d d xV= 9 l^ {±h[A+ ^R] h - |V [A + ^ 

+ hDOLTh^ + (D^)(D a h a „)} , (2.12) 
S 2 I 2 = J d d xVg~f a2 {^h 2 - V^} • 

Here h = g^h^, A = — Z) 2 , and we have freely integrated by parts. All indices are raised 
and lowered with the background metric. Furthermore, we have specified the background 
metric as the one of the <i-dimensional sphere, satisfying 

1 R 

which suffices to keep track of the expansion in the background scalar curvature terms, cf. 
eq. (2.8). All terms are generalized homogeneous in the conformal factor £ 2 and scale with 
exponents 

a = i(d-4)i/, ai = i(d-6)i/, a 2 = ±(d(l-n)-4)v, (2.14) 
for the cosmological constant, I±, and I 2 , respectively. 
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(2) 

In order to diagonalize the Hessian T k we implement the transverse-traceless (TT)- 
decomposition [56] of the fluctuations fields with respect to the spherical background, 
according to 

V = hi, + D^ v + D v ^ + D^D^a + ±g^(-D 2 a + h) , (2.15) 

for the metric fluctuations and 

= Cj + Dfi , C^tf + Drf, (2.16) 

for the ghost fields, respectively. The component fields are subject to the (differential) 
constraints 

0"*7£ = O, = 0, D^ = 0, D^ = 0, D^ = 0. (2.17) 

The resulting Jacobian determinants resulting from the TT-decomposition are exponen- 
tiated by introducing suitable auxiliary fields along the lines of the Faddeev-Popov trick. 
For this purpose, we introduce the transverse vector ghosts c T/1 , c^, the transverse vector 
bl the scalar ghosts c, c, the real scalar b, and a complex scalar s, s, which enter into the 
auxiliary action (see [13, 32] for more details). On the spherical background (2.13) it reads: 

S aux = I d*xy/7j{cZ [A - \R] c T " + ^c[A- ^R] A c 

J 1 ' (2.18) 

+ K b^ + ^biA-^Ab + SAs}. 

Substituting the decomposition (2.15), it is now straightforward to obtain the compo- 
nent field representation of (2.12) 

S% = J d d x^gl^ {^^h [A + C S R] h - [A + C T R] - d -f^R [A - \R] ^ 

- \h [(d - 1)A - R] Aa + a [^A^ _ ^R A + Aa} , 

5H 2 = I=* / d d x^g^ {^^h 2 - h^h^ - 2? [A - IR] ^ - \a [(d - 1)A - R] Aa} , 

(2.19) 

where we abbreviated 

° T= d(d-l) ' Cs= 2{d^)- ( } 
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Analogously, the quadratic form arising from the gauge-fixing term reads 

^ Sf = i / ddx ^{i^-\R^^ + i<y[{d-i)^-R] 2 ^<y] , (2.21) 

while the ghost action gives 

5 2 S ghost = V2 J d d x^t 2v {Cl [A - \R] C^ + lfj [(d - 1)A - R] A 77} . (2.22) 

Reinstalling the fc-dependent coupling constants multiplying the action monomials, it 
is then straightforward to compute the Hessian 

« r,i 1 1 5 2 T k 



(2) 



{x,y) ~ { 1)W v^v^^(*)4(v) (2,23) 

where ip = (ipi) = (h^, h, ^, a; Cj, Cj, fj, rj; cj, cj, 6j, c, c, 6, s, s) is the multiplet of all 
fluctuation fields, and (— 1)^ takes values or 1 for (pj Grassmann-even or odd, respectively. 



The matrix elements of 



F (2) 
1 k 



in field space are then summarized in the second column of 
Table 1. In order to uniformize the expressions in the gravitational sector, we introduced 
the d, n-dependent constants 

d 



c 



(d- l)(d-2) 



(d(l - n) - 2^ (1 - n) , c 2T = 2(1 - n) . (2.24) 



(2) 

Notably, the entries of T y k in the and acr sector contain the contribution from the gauge- 
fixing term only, and omit the terms originating from T™ et . The latter are subleading in a 
and can be shown to drop out of the flow equation once the Landau limit a — > is taken. 
Anticipating this result, Table 1 gives only the leading a-terms. 

The next step in obtaining the /3-functions is the construction of the matrix- valued IR- 
cutoff operator 1Z k . This operator provides a fc-dependent mass term for the fluctuation 
fields which is built from the background metric only. This implies that IZk cannot depend 
on e. The e-dependence of rj^ then enforces a modification of the cutoff schemes used 
in previous single- metric computations. Focusing on the cutoff of Type I [14] the rule for 
determining H k in a single-metric truncation adjusts TZ k in such a way that all covariant 
Laplacians are dressed by a /c-dependent mass-term according to 

A m- A + R k = P k . (2.25) 

Here R k = k 2 R^(A/k 2 ) and R^(z) is a shape function interpolating monotonously be- 
tween i?(°)(0) = 1 and lirn^oo R(°\z) = 0. For the bimetric setup of this paper, we 
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Fields 



Hessian T 



(2) 



Kernel IZk 



h T h T 
hh 



32vrG 



(d-2)(d-l) 



[^ 2ai A + £ 2oii CtR - 2A k x £ 2a ° + c 2T M^ 2q?2 ] 
£ 2Q1 A + £ 2Q1 C 5 i? - £ 2ao + c M fc ^ 2a2 



32vrGi 



(d-2)(d-l) p 



(JCJ 



s [ A - a-**]' 



(d-i) 2 

ad 2 



A 



d-1 



R 



A 



v/2£ 2 ^ [A - 



2y/2P v *=jj± 



A 



A 



I [p2 



A 2 - lRR k ] 



2V2(d-l) 
d 



p2 



V2R k 
A 2 



R p 



cc 
bb 
ss 





A- 


\R 








R k 






A- 


\R 








Rk 




d-1 
d 


A- 


dhiR 


A 


d-1 
d 


~p2 


- A 2 - 


dhRRk 


d-1 
d 


A- 


d=iR 


A 


d-1 
d 


' p2 
r k 


- A 2 - 


ahRRk 



A 



Rk 



Table 1: Matrix elements of F k ' and the coarse graining kernel to leading order in the gauge- 
fixing parameter a. The horizontal lines separate the contributions from the gravitational, ghost, 
and auxiliary sector, respectively. The explicit expression for TZ k a is provided in eq. (2.27). 



generalize this rule in the minimal sense 

A\ e=Q ^ A + R k \ e=Q , (2.26) 

i.e. Rk provides a /c-dependent mass term at zeroth order in the e-expansion. This definition 
reduces to the standard Type I cutoff implementation for the single-metric case. 

Applying (2.26) to the matrix entries Y k then determines the entries of TZk uniquely. 
The result is displayed in the third column of Table 1, with the explicit form of Vf k being 



T><T(T _ (d-i) 2 

K k - 



(Pk - ^R) 2 Pk - (A - ^R) 2 A . (2.27) 



With this result, we now have all ingredients for the explicit construction of the operator 
trace appearing on the right-hand-side of the flow equation resulting from our ansatz. 
Utilizing the block diagonal form of rj^ in field-space, this trace decomposes as 

Wfo, g] \ 9={ i + er § = s ™ + 5 o + Sgf + 5 aux . (2.28) 
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Here the left- hand- side is given by the <9 r derivative of the double expansion (2.7), while 
$2T, Sq, «Sgf, and iS aux are the operator traces constructed from the transverse-traceless 
/^-fluctuations, the metric scalar h, the gauge-fixing sector in the second block, and 
the auxiliary field contribution given in the third block of Table 1, respectively. By first 
carrying out a double expansion of the trace-arguments with respect to e, R, retaining all 
the terms indicated in (2.8), the traces can be evaluated using standard early-time heat- 
kernel techniques. Since the corresponding computation is rather technical, it has been 
relegated to Appendix B, where we also give the explicit expressions for the d- dimensional 
/3-functions. For the rest of the paper we restrict ourselves to the case d = 4 for which the 
explicit /3-functions are given in the next subsection. 

We close this subsection with a remark on the unphysical exceptional modes arising 
from working with the TT-decomposition on a spherical background. Performing a spectral 
decomposition of the component fields £ M , a in terms of A-eigenmodes, one finds that the 
two lowest scalar eigenmodes (the constant mode and the lowest non-trivial eigenfunction 
satisfying the conformal Killing equation) and the lowest vector-eigenmode (satisfying the 
Killing equation) do not contribute to and therefore require special care when eval- 
uating the operator traces (2.28). One finds that their contribution to the flow equation 
enters only at 0(R 2 ), however, so that this subtlety can safely be disregarded in the present 
computation. 

2.4 The four-dimensional /^-functions 

Based on eq. (2.28), the /3-functions 

d t gi(k) = A(g) , g = {A fc , g k , m k , X B , g B } , (2.29) 

arising in the rf-dimensional double-Einstein-Hilbert truncation are computed in Appendix 
B. To simplify our notation we will set d — 4 in the following and work with the 4- 
dimensional version of the dimensionless coupling constants (B.l) 

X k = k~ 2 At g k = k 2 G£, m k = k- 2 M k , X B k =k~ 2 A B k , g B = k 2 G B k . (2.30) 

Furthermore, we define the anomalous dimensions of the two Newton constants as 

VN = (Gt)-%G£ , V B = {Glycol . (2.31) 
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The /3-functions are most conveniently expressed in terms of the dimensionless threshold 
functions (A. 5). In this context, it turns out to be convenient to introduce a short-hand 
notation for the A- and ^-independent terms appearing in the square brackets in the first 
and second line of Table 1. The arguments of the threshold functions are then given by 
the e = 0-limit of these terms. For the h^g and /^-contributions they read 



W 2 T 



[1 - n)m k - \ k 



w 



(n - l)(2n - l)m k - \ k 



(2.32) 



Furthermore, we denote the first and second derivative of these terms with respect to e by 
w' 2T = 4n(n — 1) ink , w'q = — %n{n — \)(2n — 1) m k , 



«4r 



—An(n - l)(2n + 1) m k , 



,// _ 8 



n(n - l)(2n - l)(2n + 1) m k . 



(2.33) 



These expressions can also be obtained by setting d — 4 in (B.2) and (B.3), respectively. 

Using these notations the /3-functions in four spacetime dimensions can be summarized 
as follows. First we explicitly solve (B.ll) from Appendix B for r] N . This yields the 
anomalous dimension 



g k Bi(\ k ,m k ) 

Vn = ^-tt r • (2.34) 

The B 1 and B 2 are obtained by splitting Ft 1 ' 1 ) = (167r) _1 (_Bi + r] N B2) into the terms 
independent and linear in 1]^: 



B 



Bo 



1 = — {25 $2(^2t) + $ 2 (™o) - 5w' 2T $?(^ 2T ) - w' $1(W ) 

+ 40 w' 2T ® 3 2 (w 2T ) - 80 $l(w 2T ) - 3 <$>l + 28 <$> 3 3 + 72 <^}, 
- Y^{ 25 $1( w 2t) + $2(^0) - 5w' 2T 3>l(w 2T ) - w' <bl(w ) 
+ 40w' 2T mw 2T ) ~ 80 $*(w 2T )}. 



(2.35) 



Here the threshold functions $^(u>) and without explicit argument are understood 

to be evaluated at w — 0. Utilizing (2.34) the /9-functions for g k , \ k and m k can be obtained 
from equation (B.13): 



Pg = [ 2 + Vn ] 9 k , 



X =(Vn ~ 2)\ k + -g k [4F( 2 >°) + 2(2n - 1)F™] , 

71/ 

2tt 



Pm =(Vn - 2)m k + 



n(l — n) 
14 



g k [2F^ - F™] . 



(2.36a) 
(2.36b) 

(2.36c) 



The running of the background couplings is governed by 
f3 gB = [2 + r%] g% , 



(2.37a) 



/3 AB = fo* - 2)A* + [8(1 - n)F^ + AF^ - 2(3 - 2n)F^] , (2.37b) 

where the background anomalous dimension is given by 

V * = I6n (F^ - F^) . (2.38) 
The general result for the functions F^' j \g k , \ k , m k ) is given in Appendix B. Setting d = 4 



it simplifies to: 
^(o,o) ; 



^(1,0) 



^(2,0) 



(4tt) 2 
1 

: (4^j2 

1 

: (4tt) 2 
1 

: (4^)2 



[5g 2 1 (w 2T ) + g 2 1 (u'o)-4$^ , 



|gi(w 2T ) - f q 2 2 (w 2T ) 



i „i 



12 2 3 ^1 



W 3 (w 2T ) + 2q 2 3 (w ) - 5 w' 2T q 2 2 (w 2T ) - w' q 2 2 (w ) + 24 <^ + 12 $ 2 j , 
- | qf(w ) - | w 2T (w 2T ) + hq 2 (w ) + fqi(w 2T ) + f w' 2T ql{w 2T ) 



~ ?9 3 S («*r) 



1 $2 + U $ 3 + 12$ 4 



'(4tt) ; 



30g|(w 2 T) + 6g|(w ) - 10gf(u> 2 T) - 2g 2 (u? ) - 20 w' 2T ql(w 2T ) - 4w' ql(w ) 

+ 5 (w' 2T ) 2 q 3 2 (w 2T ) + K) 2 q 3 2 (w ) - §< T q 2 2 (w 2T ) - ±< g 2 > ) - 36 $^ - 480^ 

(2.39) 

The /9-functions (2.36) and (2.37) together with (2.34) and (2.38) constitute the main result 
of this section. Their properties, in particular the fixed point structure they give rise to, 
will be investigated in the next section. 



2.5 Bimetric vs. single- metric truncations 

At this point it is illustrative to compare the bimetric ansatz (2.2) to the single-metric 
(SM) Einstein-Hilbert truncations studied previously [8, 14, 20, 21, 25, 31] where 

1 



16ttGI m 



d d x^~g [R - 2AH 



(2.40) 



It is this functional which, for the Einstein-Hilbert case, corresponds to T k [g] in eq. (2.2). 
It depends on the background only via g^ u = g^ u + h^ u , that is, it has no extra background 
dependence: T^ ct [g,g\ = T k [g]. 
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How is the functional (2.40) and its predictions for the RG flow related to those of the 
bimetric ansatz (2.2)? The correct mapping of the bimetric computation onto its single 
metric analog consists in simply omitting the second and third integral on the RHS of 
(2.2). This amounts to setting 1/Gf = 0, Aj? = 0, M k = and identifying 

Gt M = Gt Af 1 = At (2.41) 

While this mapping scheme is very natural, the reader might argue that there is another 
one that seems equally plausible. In a certain sense, the single metric truncation does not 
distinguish g^ u from g^ u , so one could be motivated to keep all three integrands in (2.2), 
but replace g^ u with the expectation value metric g^ u everywhere. The result is a single 
metric functional of the form (2.40) with the couplings 

1 1 1 Af 1 _ ^ At M k 

+ i nSM ~~ + ~Fm ~ nA ■ 1^-42) 



Ul Or,, Or, Or, Or, Or, Or, 

We must emphasize that the identification (2.42) is actually not the correct way of relating 
the new bimetric calculation to their old single metric counterpart. The reason is that 
if one proceeds in this way one has to retain the second and third integral of (2.2), with 
g = g, also on the RHS of the flow equation. The identifications (2.42) amount to setting 
g = g directly in the ansatz, i.e. prior to computing the Hessian. As a result, T k 2 \ in this 
case, receives contributions also from the second and third integral of (2.2), involving the 
background couplings G7j?, Af , and M k . Since those contributions are absent in the actual 
bimetric calculation it is clear that (2.42) cannot be correct. 

Within our present computational setting the /3-functions for the dimensionless single- 
metric couplings {g SM , A SM } are easily derived noticing that the RHS of the single-metric 
RG flow equation is given by the zeroth order terms F^°'°\g SM , A SM ) and F^^g™, A SM ) 
in the e-expansion of (B.7). Thus, the equations read 

(3 gSM = [2 + ^}gl M , (2.43a) 
/3 a sm = fa™ - 2)A* M + 8ngl M F^ , (2.43b) 

where ?7^ M = lGirgf^F^ ' 1 ^ is the anomalous dimension of G7| M . We shall come back to 
these /3-functions shortly. 
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3 RG flow of the double-Einstein-Hilbert truncation 



We now investigate the RG flow resulting from the /3-functions (2.36) and (2.37). The 
system without mixed term, m k = 0, is analyzed in Subsection 3.1 while the properties of 
the full system are discussed in Subsection 3.2. 

3.1 The double-Einstein-Hilbert truncation without mixed term 

In order to understand the properties of the bimetric RG flow, we first discuss the simpler 
four-parameter truncation without the mixed term oc (vWy 7 ^)™- ^he corresponding j3- 
functions can be recovered from eqs. (2.36) and (2.37) which requires some care though. 
(See the remarks at the end of Appendix B.) In d — 4, the system of /3-functions reads 

(3 g = [2 + m }g k , (3.1a) 

f3 x =( VN -2)\ k + 4ng k F^°\ (3.1b) 

^ = [2 + 1^] gl (3.1c) 

/3 AB =(V B N - 2)Af + 8irg* [F^ - \ F™] , (3.1d) 

where tjn and i]® are the anomalous dimensions defined in (2.34) and (2.38), respectively. 
This system is decoupled in the following sense. Eqs. (3.1a) and (3.1b) close among them- 
selves and are independent of the background couplings g B and A B . Once a solution of this 
subsystem is given, we can insert it into the remaining /3- functions (3.1c) and (3. Id) and 
solve the resulting two differential equations for g B and A B . To find fixed point solutions in 
this four-parameter system, it is enough to search for fixed points ((?*, A*) of the subsystem 
(3.1a) and (3.1b), substitute (g*, A*) into (3.1c) and (3. Id), and then look for zeros (gf, A B ) 
of the background /3-functions. 

Following this strategy and using the optimized shape function [57], a numerical search 
unveils the following fixed point structure. The subsystem (3.1a)-(3.1b) admits a Gaussian 
fixed point (GFP) at A*) = and a NGFP with A*) ^ 0. Each of these fixed points 
gives rise to a pair of zeros of the background /3-functions (3.1c)-(3.1d). One of them 
corresponds to a background GFP with (g B , A B ) =0, the other to a background NGFP at 
(gf, A B ) 0. The resulting four combinations of fixed points are summarized in Table 2. 

Having found a fixed point g* of the four-parameter system g = (gj) = (g, A, g B , A B ), we 
compute the stability matrix = <9j/3j(g*) which governs the RG flow linearized around 
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the fixed point: 

kd k gi(k) = Bij (gj(k) - . (3.2) 
Setting t = \n(k), the general solution of (3.2) reads, in the nondegenerate case, 

Sj(k) = g*, + r « e " n e ~ 6nt V i • ( 3 - 3 ) 

n 

Here V n are the right eigenvectors of the stability matrix, with eigenvalues — 9 n , and 
r n e ian = C n are free constants of integration. They can be complex except when 9 n 
happens to be real (then a n = 0). Critical exponents with Ke(9 n ) > correspond to 
relevant scaling fields. They grow when k is lowered, i.e. they amount to UV attractive 
directions. 

Based on the results shown in Table 2, the stability matrix governing the linearized 
flow near the NG-NG-FP can be evaluated numerically. Upon its diagonalization we find 
a pair of complex conjugate critical exponents Q\ = 9* 2 = & + iO" with 

& = 4.468, 9" = 4.240, (3.4) 

together with the background critical exponents 

#4 = 4 and #5 = 2. (3.5) 

Since all critical exponents are positive, the NG-NG-FP is UV-attractive in all directions. 
While 9' and 9" depend on the cutoff chosen, our results for 6> 4 and 6> 5 are universal, i.e. 
cutoff independent. This follows from the special structure of the system (3.1), which 
implies a stability matrix which has a lower triangular form. Consequently, the critical 
exponents of the background /3-functions (3.1c) and (3. Id) are given by: 



6 a = -^- = 2-t%, and 6l 6 = -^ = -2(1 + ^). (3.6) 



At a non-Gaussian fixed point the background anomalous dimension function is 77^ = —2, 
establishing that the critical exponents (3.5) are indeed universal. 

At this stage, it is useful to pause and have a closer look at the general mechanism that 
generates the NGFP for the background couplings. For a typical background coupling 
■uf f d A x*JgO with mass-dimension d m the /3-function for the corresponding dimensionless 
coupling = k~ dm u® will assume the form 

dt u B k =-d m u B k +f(u£,uD, (3.7) 
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Fixed Point 


(?* A* 


9? ^ 




G-G-FP 
G-NG-FP 
NG-G-FP 
NG-NG-FP 




1.055 0.222 
1.055 0.222 



2.20 -0.131 

-41.649 0.578 



-0.29 
0.234 
0.234 -24.06 



Table 2: This table shows all fixed points occurring in the RG flow of the four-parameter system 
without mixed term given by (3.1). The products g* A* and gf\® are also given. 

where and «™ lx denote the dimensionless coupling constants multiplying the interaction 
terms including the "genuine" metric and a mixture of g and g, respectively. The later are 
determined through their corresponding /3-functions so that /(m^,m™ 1x ) is a fixed number. 
The fixed point value for the background coupling is then obtained by solving the linear 
equation (3 u b = 0: 

u*=(d m )- 1 Hu$ > u? x ). (3.8) 

Notably, this mechanism works for all dimensionful background coupling constants, but 
fails if is power-counting marginal, i.e. it has d m = 0. In the latter case does not 
obtain a finite fixed point value and runs logarithmically in the UV. Clearly, it would be 
highly desirable to get a better understanding of this very perculiar feature. A complete 
clarification of this issue is, however, beyond the scope of the present work and will be left 
for future study. 

A further question arising naturally at this point is how the bimetric results relate to 
those coming from the single-metric truncation. We must then compare the subsystem 
(3.1a)-(3.1b) to the system of /3-functions arising in a single-metric truncation as given in 
(2.43). We find that the system (2.43) has a NGFP located at 

g* M = 1.129, Af 1 = 0.216. (3.9) 

Diagonalizing the stability matrix yields to complex conjugate critical exponents with 

^ M = 1.709, 0g M =3.44. (3.10) 

These figures refer to exactly the same computational setting (cutoff type, shape function, 
field parametrization, etc.) as the bimetric computation above. According to the discussion 
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(a) Bimetric truncation 



(b) Single-metric truncation 



Figure 1: Phase portrait of the bimetric (left) and single metric truncation (right). The phase 
portraits show that both systems are qualitatively very similar. The plots also show the boundary 
of the physical parameter space. 

that led to the identification (2.41) we expect that, if the single metric truncation is a 
reliable approximation to the bimetric truncation, we should find gf M ~ g* and A^ M ~ A*. 
Comparing the single metric NGFP with the NG-G-FP or the NG-NG-FP of Table 2 we 
see that these relations are satisfied remarkably well. However, the critical exponents of 
the two systems are rather different. This confirms that the bimetric corrections are indeed 
important at the quantitative level. 

We close this subsection by comparing the phase portraits resulting from solving the 
bimetric (3.1) and single-metric flow equations (2.43). An illustrative set of numerically 
obtained sample trajectories is shown in Figure 1. Remarkably, both flows exhibit the 
same qualitative behavior, despite being based on two quite different systems of differential 
equations. In both cases we can distinguish trajectories that run for decreasing k towards 
negative cosmological constants, referred to as trajectories of Type la in the terminology 
of [21], to positive cosmological constant (Type Ilia), and to a vanishing A. The latter, 
single trajectory (of Type Ha) is a separatrix; it crosses-over from the NGFP to the GFP. 
Both in the single metric and the bimetric case all trajectories in the upper half plane 
(positive Newton constant) are pulled into the non-Gaussian fixed point as we send k — > oo. 
The shaded region in the plots of Figure 1 is delimitated by a line where the anomalous 
dimension t]n diverges [21] and does not belong to the physical parameter space. In the 
bimetric case, too, all trajectories of Type Ilia terminate at a finite value of k at this 
boundary. Exactly as in the single-metric case they cannot be continued to the physical 
point k = within the truncation used [26]. 
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Fixed Point 




9* A * 


9*K 9?^ 


n = 2 


UV-NGFP 


1.273 0.237 0.025 


2.009 0.133 


0.3017 0.267 


IR-NGFP 


0.821 0.210 -0.080 


-1.339 0.251 


0.172 -0.336 



Table 3: Position of the fixed points occurring in the double-Einstein-Hilbert truncation (2.2) 
for the sample value n = 2. The products g* A* and g®\® are also given. 

3.2 The double-Einstein-Hilbert truncation including mixed in- 
teractions 

Now we switch on the mixed interaction term oc and search for zeros of the complete 
system of /3-functions, that is, for fixed points g* = A*,m*, , \®). We analyze (2.36) 
and (2.37) for general values of the parameter n = 2,3,4. .. . Notably, this system has 
the same lower triangular structure as in the case described in Subsection 3.1. Thus 
the critical exponents corresponding to the background couplings are exactly the same. 
However, allowing for a non-zero value of entails that the NGFP found before now 
splits into two different NGFPs, referred to as UV-NGFP and IR-NGFP, respectively. 

We now exemplify these structures for the special case n = 2. The corresponding 
position of the two fixed points is given in Table 3. Linearizing the /3-functions near the 
fixed points, we find that the flow near the UV-NGFP (IR-NGFP) is governed by a pair of 
complex conjugate critical exponents #i = Q\ = ®' + W ■> a positive (negative) real critical 
exponent # 3 , and the same universal background critical exponents as before 

UV-NGFP: #' = 4.800 #" = 8.722, # 3 = 8.151, # 4 = 4, 6 5 = 2, (3.11a) 
IR-NGFP: #' = 3.131 #" = 2.821, # 3 = -24.814, # 4 = 4, # 5 = 2. (3.11b) 

These stability properties then motivate the denominations "UV-NGFP" for the fixed 
point whose eigendirections are all UV-attractive and "IR-NGFP" for the one with one 
UV-repulsive eigendirection, respectively. In fact, we shall shortly see, that there exists 
a cross-over trajectory which emanates from the UV-NGFP in the UV and ends at the 
IR-NGFP in the limit k -> 0. 

The position and critical exponents of the UV-NGFP and IR-NGFP for general values 
n are shown in Figure 2. Notably, both the UV-NGFP and IR-NGFP occur for all values 
n > 2. Interestingly, the fixed point value m* tends rapidly to zero as we increase the value 
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5 10 15 20 5 10 15 20 5 10 15 20 

n n n 

(a) UV- and IR-g* (n) (t>) UV- and IR-A, (n) (c) UV- and IR-m* (n) 




(d) UV- and IRS' (n) (e) UV- and IR-6" (n) (f) UV- and IR-fl 3 (n) 

Figure 2: The n-dependence of the two NGFP in the {g, A, raj-subsystem (a)-(c), together 
with their n-dependent stability coefficients (d)-(f). The plotted circles represent the IR-NGFP 
whereas the "*" represent the UV-NGFP. 

of n, in both cases (cf. Figure 2c). In Figure 3a we show the 3-dimensional phase diagram 
of the {g, A, m}-subsystem with its UV-NGFP and IR-NGFP. The sample trajectories 
in this plot have been obtained numerically. Figure 3b depicts the projection of the 3- 
dimensional flow onto the g-X— plane. Remarkably, this projection looks very similar to 
the flow generated by the single metric Einstein-Hilbert truncation shown in Figure lb. The 
Figure 3a allows us to identify the 3-dimensional generalizations of the familiar trajectories 
of Type la and Ilia, i.e. their "lift" to g-X-m— space. The g- X— projection of the 3- 
dimensional trajectories displayed is almost identical to what one obtains from the single 
metric truncation. In particular, according to both the 3-dimensional and 2-dimensional 
system of RG equations the trajectories of Type Ilia terminate at finite k in a boundary 
singularity. For both the 3-dimensional and the 2-dimensional single-metric RG equations 
this class contains trajectories whose turning point (the point where (3\ = 0) is arbitrarily 
close to the GFP. Those trajectories spend a very long RG time near the GFP. As a result, 
their termination at the boundary singularity can be deferred until very late RG times, that 
is, small scale k. All these features common to the 2- and 3-dimensional RG flow fortify the 
rather impressive robustness of the projected phase portrait. In refs. [58] and [51] it has 



22 



(a) Three-dimensional phase portrait 
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2.5 
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0.05 0.10 0.15 0.20 0.25 0.30 0.35 
X 

(b) Two-dimensional projection 

Figure 3: (a) Phase portrait of the {g, A, m}-subsystem. (b) The projection onto the g-X— plane 
of some sample trajectories. The special trajectory crossing over from the UV-NGFP to the 
IR-NGFP is shown in both plots. 

been pointed out that if one uses a trajectory of Type Ilia to define the quantum theory, 
the mechanism of "using up" a long RG time near the GFP is instrumental in obtaining 
a theory with a classical regime given by General Relativity where Gj~ = constant for a 
long interval of scales k. The same argument also applies to the "lifted" trajectories on 
g-X-m— space. For a tentative matching of these model trajectories to the observations 
made in Nature the arguments in [58] and [51] basically remain unaltered therefore. 

A highly intriguing new feature of the 3-dimensional flow is the emergence of a novel 
IR-NGFP and, as a result, the cross-over trajectory connecting the UV-NGFP to the 
IR-NGFP. It emanates from the UV-NGFP in the UV and for k — > approaches the IR- 
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NGFP along its only IR attractive direction, see Figures 3a and 3b. For no value of k this 
trajectory gets close to the GFP. As a result, it does not give rise to a classical regime: 
the dimensionful parameters A£, M k have a significant scale dependence all the way 
down from u k = oo" to k — 0. 

While the cross-over trajectory in the form obtained within the double-Einstein-Hilbert 
truncation cannot be used as a model for the real world, its existence and, more generally, 
the emergence of an IR fixed point, is a highly welcome feature of the new truncation. 
One reason is that, for the first time, we have found a trajectory which has positive A 
everywhere and which does not terminate in a singularity. As to yet it had always been 
unclear what kind of average action could avoid the singularity near A = 1/2. (See [26] for 
an early attempt at solving this problem.) Moreover an IR fixed point, albeit at a different 
location, had been argued to lead to a particularly realistic late-time cosmology and an 
explanation of the "recent" cosmic acceleration [52]. 

Clearly, for phenomenological purposes the IR-NGFP is not yet satisfactory. But the 
fact that it occurs at all in the new truncation is very encouraging. It shows that we are 
"on the right track" towards a better understanding of the infrared. 

4 Discussion and Conclusions 

In this paper we analyzed the functional RG equation for full-fledged Quantum Einstein 
Gravity within the framework of a novel class of approximations, the bimetric truncations 
of theory space. In the ansatz for the effective average action we included two different 
Einstein-Hilbert actions, one for the dynamical and another for the background metric, 
as well as a non-derivative term which mixes the two. Our main interest was in assessing 
whether the resulting RG flow lends itself to the Asymptotic Safety construction of a 
microscopic theory of quantized gravity. One key result is that as far as one can tell 
within the restrictions of the approximation, the answer is clearly positive. In particular, 
we discovered a bimetric generalization of the non-Gaussian fixed point known to exist 
within the single-metric Einstein-Hilbert truncation. This confirms earlier expectations 
[43] originating from a similar bimetric computation within conformally reduced gravity. 
Moreover, we also found a second non-trivial fixed point which might control the theory's 
IR behavior. 
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The main technical innovation of the present paper is the conformal projection tech- 
nique. As its main virtue it allows to probe (a part of) the bimetric gravitational theory 
space without the necessity of evaluating complicated operator traces involving the two 
Laplacians constructed from both the "genuine" metric g and the background metric g, 
Tr[/(— D 2 , — D 2 )}. A complementary approach could be based on expanding Tk[h;g] in 
terms of the n-point functions of h^ u = g^ u — g^ u . Results for full-fledged quantum gravity 
using the latter approach will be reported elsewhere [59]. 

It is impressive to compare the /3-functions for the couplings g, A, (3.1a) and (3.1b) to the 
corresponding ones arising in the single metric computation, eqs. (2.43). The latter use the 
same gauge fixing and auxiliary field construction, so that this can serve as an illustration 
which makes the new features of the bimetric computation transparent. Obviously, the 
/9-functions (3.1) and the single metric result (2.43) are very different in their analytical 
structure. In particular all coefficients multiplying threshold functions, as well as their pole 
structure at A = 1/2 and A = 3/4, are manifestly different. In this light it is miraculous that 
the single-metric and bimetric treatments, nevertheless, give rise to a non-Gaussian fixed 
point with very similar properties. For this reason we may be optimistic that the bimetric 
generalizations of the RG flows we know already will, at least qualitatively, confirm the 
essential features of the single-metric truncations. This concerns in particular the viability 
of the Asymptotic Safety program. 
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A Threshold functions 

In this appendix, we collect various definitions and review the central properties of the 
threshold functions used in the main text. In order to evaluate the trace contributions, we 
use the early time expansion of the heat kernel. For the operators appearing in Table 1 



25 



this can be done utilizing the the "master formula" 



Tr s [W(-D 2 )} = ^tr s (l s ) ^Q d/2 [W] J d d x^§ + ^Q d/2 -i[W] J d d x^R 



+ 



(A.l) 



where 



1 f°° 

Q n [W} = — - / dzz^Wiz), (A.2) 
1 \ n ) Jo 

and the higher-order terms outside the truncation are indicated by the dots. The subscript 
s = 2T, IT, indicates that the Laplacians act on transverse-traceless symmetric tensors, 
transverse vectors, and scalars, respectively. The algebraic trace tr s (l s ) counts the degrees 
of freedom in the corresponding sector. Setting d s = tr s (l s ) we have 



d 2T = i(d-2)(d+l), d lT = d-l, 



do 



(A.3) 



When performing the double-expansion of the traces (2.28) in e, R, the Q-functionals ap- 
pearing in the trace-evaluation can be related to the standard threshold functions via 

&dt(Z k R k ) 



Qn 



Z k {z + R k + wk 2 y 



T(n+q) , 2(n+g-p+l) 
r(n) 



2& n+q (w)-ri N $ p n+q (w) 



, n + q > . 

(A.4) 

Here r/ N = —d t \nZ k and z = — D 2 . The cutoff-scheme dependence is encoded in the 
standard dimensionless threshold functions [8] 



1 



dzz 



n _ l R^(z)-zR^'(z) 



r(n) 



dzz 



n-l 



z + R(°\z)+w 



p ' 



(A.5) 



+ R(°\z) +w 



p ■ 



To simplify the notation of the traces including the running Newton's constant Z k) the 
right-hand-side of (A.4) suggests introducing 



(A.6) 



As we shall see later on, this short-hand notation will considerably simplify the notational 
complexity of the evaluated traces. 

For the explicit evaluation of the flow equations, we resort to the optimized cutoff [57]. 
In this case the shape-function entering into (A.5) is given by 



Rl°)(z) = (1-3)0(1-*) 



(A.7) 
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and the integrals appearing in the threshold functions can be carried out analytically: 

* iW = r(,'+i)(n-'»)> ' ^''- r^Hi+W (A ' 8) 

We exclusively resort to this type of cutoff in the main part of the paper. 



B Double-Einstein-Hilbert truncation: deriving the 
^-functions 

In this appendix we present the technical details entering into the evaluation of the trun- 
cated flow equation (2.28). We start with the explicit calculation of the operator traces 
constructed from the entries of Table 1 using the intermediate formulas reviewed in Ap- 
pendix A. Subsequently, the projection of the RG flow gives rise to our main result, the 
rf-dimensional /3-functions of the double-Einstein Hilbert truncation, (B.13) and (B.14) 
below. 

Our starting point are the operator traces entering into (2.28). Here we first focus 
on the gravitational sector sourced by the transverse-traceless h^ v and the trace h, before 
evaluating the universal contributions from the gauge- and auxiliary sectors given by the 
second and third block of Table 1, respectively. All expressions are given in terms of the 
dimensionless coupling constants 

\ k = k- 2 A£, g k = k d - 2 G%, m k = k~ 2 M k , \ B k =k' 2 A B k , g B = k d ~ 2 G B . (B.l) 

Comparing the first and second line of Table 1, we observe that the expressions enter- 
ing into S and S 2 t are structurally very similar. This motivates the following helicity- 
dependent definitions, in terms of which the evaluated traces take the same structural 
form. First, setting £ = 1, we identify the arguments that will enter into the dimensionless 
threshold functions as the A- and ^-independent parts of the square brackets in Table 1, 

w 2 t = c 2T m k - 2A fc , w = c m k - ^\ k . (B.2) 

The (d, n)-dependent constants c 2 t and c given in (2.24). Furthermore, the e-expansion 
motivates introducing short-hand expressions for the e-derivatives of the m, A-dependent 
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terms in Table 1 

W ' 2T = c 2 t « 2 m k - 2 a X k , u>g = c a 2 - ^ty aoAfc , 

w 2T = c 2 Ta 2 (a 2 - l)m k - 2a («o - l)A fc , = c a 2 (tt2 - 1)m k - ^a (a - \)\ k , 

(B.3) 

with ctj defined in (2.14). In terms of these <So and S 2 t assume the same, though spin- 
dependent, form. The explicit computation, utilizing the formulas in Appendix A yields 

S s = J ^ m q 1 d/2 (w s )k d J d d x^g 

+ [hd/2-i( w s)-C s q 2 d/2 (w s )]k d - 2 J d d x^R 

k d [ d d x^~§e 



d s 
(4n) d / 2 



-V^d/2+iK) + w s q d/2 {w s ) 
^(^rr + C s )q% 2 {w s ) + \w' s q% 2 _ 1 {w s ) 



(B.4) 



- f [vd{d - Q)q% 2+1 {w s ) +Aw' s q 3 d/2 (w s )) ] k d ~ 2 J d d x^R 



+ 



(47r) d / 2 



v 2 d(d+2) , 



id - Q) 2 q% 2+2 {^s) + W-GXti/z+iM + K)^/ 2 (™ 



vd 



16 



v(d-6)d 
16 



((d - 6)1/ - 2) gj /2+1 K) - |< gj /2 K) A; d / d d x^ge 



The gauge-sector, given by the second block in Table 1, likewise leads to 



+ 



+ 



(4n) d / 2 

vd 
(4n) d / 2 

v 

(4n) d / 2 



^/2^y ^VI + Wd/2-ik d - 2 J d d x^R] 
(d + 2)<t>* /2+2 + (d - l)<^ /2+1 ] k d J d d x^e 

(d-l)(d? -2d- 12) ^2 , d 3 +9rf 2 -34rf+12 ^s3 , 3^+2)^4 

6d ^d/2 "T 6 (d-l) ^d/2+1 "T 2 (d-l) ^rf/2+2 



(B.5) 



A; d - 2 / d d x^Re 



~ 2(£w [(1 - ")(<* " 1)^/2+1 - (1 - 21/ + d*/)(d + 2)$ 
+ §i/(d + 6)(d + 4)(d + 2)$^ /2+4 ] A; d | d^e 2 . 



,3 

d/2+2 



Here all the threshold functions are evaluated at w — 0, and we dropped the argument for 
notational simplicity. 

The auxiliary fields contribute to the running of the background couplings only. Their 
trace is universal, i.e. it does not depend on any coupling constant. It reads 



^k d J d d xV§[(d-l)^ d/2 + (^ /2 _i 



+ d(d-l) ®d/2 ) K K 



(B.6) 
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Also here, all threshold functions are evaluated at zero argument, which is dropped for 
simplicity. 

When constructing the /3-functions from the ansatz (2.1), it turns out to be useful to 
organize the right-hand- side of the flow equation in terms of the following double expansion 
in R, e: 

dMg,g}\ 9 = {1+ er s = E F m (9k,\k,m k )k d ! d d x^g e k {Rk~ 2 ) 1 . (B.7) 

k,i 

The dimensionless expansion coefficients F^ k ^ depend parametrically on n, d as well as on 
the couplings g k , \k and m k . They are, however, independent of the background coupling 
constants <?f,A^. Their explicit form is easily read off by substituting (B.4), (B.5), and 
(B.6) into (2.28) and comparing to the expansion (2.28). We find 



F(0 ' 0) = (4^72 [ rf2T 9d/2(^2T) + q\,*(wo) - d& d/2 ] , 



^(0,1) 



( 47r )d/2 



rf 2T {5^/2-1(^21?) - C T q d/2 (w 2T )} + |?d/2-i(^o) - C s q 2 d/2 (w ) 



d 2 -d+l ^,2 _ dfol 



( 47r )d/2 

_ ud(d-6) 
4 



- d 2T w' 2T q d/2 (w 2T ) - w' q 2 d/2 (w ) 

{ d2T9d/2+i ( w 2t) + q 2 d/2+ i(wo)} + i/d(d + 2)$3 /2+2 + I/d(d - 1)^/2+1 



(B.8) 



together with 



( 47r )d/2 



- ^ {4t(^ + CV^^t) + (^ + C 5 )^ /2 K)} 
- | {d 2 T w' 2T ^ /2 _i(w 2T ) + qj/^iwo)} + 2d 2T C T w 2T q^ /2 (w 2T ) 

+ 2C S w' q d/2 {w ) + {4^^/2+1 (^ 2T ) + Cs^ /2+1 (w )} 

K^-l)(rf 2 -2d-12) ^2 , -, ^+9^-34^+12 ^3 , 3d(d+2) ^ 

H 6d V ^/2 + ^ 6(5=1) V rf/2+l + ^ 2(5=1) V d/2+2 



(B.9) 



^(2,0) 



( 47r )rf/2 



" 2Ci( 16 +2) (d - 6) 2 {rf2T^/2+2(^2T) + <?d/2+ 2 K))} 



+ 



i/d(d-6) 



{rf 2T w' 2T q d/2+1 (w 2T ) + w' q d/2+1 (w ) } 



+ d 2T {w' 2T ) 2 q 3 d/2 {w 2T ) + {w' ) 2 q 3 d/2 {w ) - \ {d 2T w 2T q d/2 {w 2T ) + w'q q% 2 {w G )} 

■((d - 6)U - 2) { 6^/2+1 (^2t) + ^/ 2 +l(™o)} 



ud(d-6) , 
16 

1)^.2 ^(l-2i/+di/)(d+2)^3 ^ 2 (c(+6)(c(+4)(d+2) ^5 
2 ^d/2+1 2 V d/2+2 4 ^<i/2+4 



(B.10) 
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Note the implicit r/Ar-dependence contained in q^iw). 

Substituting (2.7) into (B.7) we obtain a system of coupled differential equations gov- 
erning the RG dependence of the dimensionless couplings (B.l). The anomalous dimension 
for Newton's constant is found solving the linear equation 

VN = lQng k F^ (B.ll) 

with F^ 1 ' 1 ) given by (B.9), for r] N . The result has the structure 

_ fffc^i(A fc ,m fc ) 
Vn - z 5-77 r , (B.12J 

1 - 9k -D2lA fc ,/i fc J 

where the functions Bi and B 2 can be read off from (B.9) by splitting the expression for 
F( 1 ' 1 ) into the terms independent of, and linear in t)n (as in [8]). 

Based on this solution of r/ N , the j3- functions for the dimensionless g k , \ k and m k close 
among themselves: 

/3 g = [ d - 2 + rj N } g k , 

fix =( VN - 2)\ k + 8ng k ^ [2(d - 2)F™ - (2 - nd)F^} , (B.13) 
fi m =( VN - 2)m k + 16ng k [(d - 2)F^ - F^} . 

The /^-functions governing the running of the background couplings are 
fo = [d-2 + r&] gf, 

/3 AB = ( v * - 2)A* + 8ngf ^ \^f F im + 2(d - 2)F( 2 > ) - (d + 2 - nd)i^ . 

(B.14) 

The background anomalous dimension id linear in gf and given by 

rj* = 16n (FM - F™) gf . (B.15) 

This completes the derivation of the /3-functions of the bimetric Einstein-Hilbert trunca- 
tion including the non-trivial mixed term proportional to m k . They are valid for arbitrary 
space-time dimension d, and constitute the central result of this appendix. Restricting to 
d — 4, their properties are discussed in Section 3. 

We close this appendix by recovering the /3-functions for the double-Einstein-Hilbert 
truncation without mixed terms oc M k , eq. (3.1). These readily follow from the general 
result (B.13) and (B.14) by setting m k = and eliminating i^ 2,0 ) via the constraint 

(d-2)^ 2 '°)-F< 1 '°)| mfc ^ = 0, (B.16) 
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which results from the consistency condition /3 m | m=0 = 0. Restricting to d — 4, one then 
arrives at the /3-functions (3.1) analyzed in Subsection 3.1. 
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